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. Abstract. We construct identities of Pohozhaev type, in the context of elastostatics and 
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ir~> ■ 1 Introduction 

Identities of Pohozhaev type have been widely used in the theory of partial differential equations, 
in particular for establishing non-existence results for large classes of forced elliptic boundary 
q \ value problems and eigenvalue problems, [201 EH G2]. The purpose of this note is to obtain and 

apply analogous identities in elastostatics and elastodynamics, which have not (to the authors' 
knowledge) been developed to date. Our approach will be based on a fundamental identity 
first introduced by Noether in her seminal paper |16j that connected symmetries of variational 
^ ' problems to conservation laws of their Euler-Lagrange equations. 

As noted in [19] . the identities originally due to Pohozhaev, [201 [21], owe their existence to 
Noether's identity. For classical solutions of the linear equation Au+Au = such an identity was 
obtained by Rellich in |24] . Further, in [25] . Rellich established an integral identity for a function 
belonging to certain function spaces, without any reference to differential equations it may 
satisfy. The Rellich identity has been generalized by Mitidieri, [121 113] , for a pair of functions. 
General Rellich-type identities on Riemannian manifolds have been recently established in [51 [6] 
by use of Noether's identity applied to conformal Killing vector fields. 

In [20] . Pohozhaev established an integral identity for solutions of the Dirichlet problem 
for the semilinear Poisson equation Au + A/ (it) = in a bounded domain with homogeneous 
Dirichlet boundary condition. Later, for solutions of general Dirichlet problems, he obtained 
in [21] what is now called the Pohozhaev identity. Such identities became very popular after 
the paper of Pucci and Serrin, [22], where, on p. 683, the relation with the general Noetherian 
theory is mentioned. See also the earlier paper by Knops and Stuart, jllj . and remarks in the 
second author's 1986 book [19]. The Noetherian approach to Pohozhaev's identities was further 
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developed and applied in [3j HI [231 EH • Additional applications of Rellich-Pohozhaev estimates 
to nonlinear elliptic theory can be found in [26} I27j . while applications to nonlocal problems 
appear in [8]. With regard to geometric applications, [H [91 [28] develop a systematic approach 
to Pohozhaev-type obstructions for partial differential equations invariant under the action of 
a conformal group. For a relation between the Lie point symmetries of the nonlinear Poisson 
equation on a (pseudo-) Riemannian manifold and its isometry and conformal groups see [2]. 

In dynamical problems, the conformal invariance of the wave and Klein-Gordon equations 
was used by Morawetz, [14] . to establish several very useful integral identities. These were 
applied by her and Strauss, [291 [30], to the study of the decay, stability, and scattering of waves 
in nonlinear media. In the final section, we will generalize Morawetz' conformal identity to some 
dynamical systems governing waves in elastic media. Applications of our identity to decay and 
scattering of elastic waves will be treated elsewhere. 

In elastostatics, the independent variables x € M. n , for n > 2, represent reference body coor- 
dinates, while the dependent variables u = u(x) = . . . , u n (x)) represent the deformation 
of the point x. The independent variable x will belong to a bounded or unbounded domain 
C R n that has sufficiently (piecewise) smooth boundary d£l. We use v to denote the outward 
unit normal on d^l. For elastodynamics, we append an additional independent variable, t, repre- 
senting the time, and so u = u(t,x). The partial derivatives of a smooth (vector) function u(x) 
are denoted by subscripts: 

k _ du k k _ du k k d 2 u k 

Ui : " ~dx~' Ut : " ~dt* U v - dx~dx~> etC ' 

The n x n spatial Jacobian matrix Vu = (u k ) is known as the deformation gradient. 

We shall consistently use the Einstein summation convention over repeated indices, which 
always run from 1 to n. We assume that all considered functions, vector fields, tensors, functio- 
nals, etc. are sufficiently smooth in order that all the derivatives we write exist in the classical 
sense. When we say that a function is "arbitrary", we mean that it is a sufficiently smooth 
function of its arguments defined on the domain 0. Extensions of our results to more general 
solutions will then proceed on a case by case basis. 

2 Noether's identity 

A vector field 

v = ?(*,«) A + ^ (a . u )^_ 

on the space of independent and dependent variables induces a flow that can be interpreted as 
a (local) one-parameter group of transformations. The vector field is known as the infinitesimal 
generator of the flow, [J2] . For example, the particular vector field 

9 , i 9 
v = axi- h bu -—, 

OXi ou % 

where a, b are constant, generates the group of scaling transformations 
(x,u) i—>. (X a x,X b u). 

The action of the group on functions u = f(x) by transforming their graphs induces an action 
on their derivatives. The corresponding infinitesimal generator of the prolonged group action 
has the form 

pr«v = C (x, u)-£-+ f(x, U ) A + $(s, u, Vn)^-, (1) 



Pohozhaev and Morawetz Identities in Elastostatics and Elastodynamics 



3 



where 



+ 



dxj du k 3 



dx. 



9C k i i 



(2) 



and Dj = d/dxj + Ujd/du k denotes the total derivative with respect to Xj. See [19] for a proof 
of this formula, along with its extension to higher order derivatives. 
For a first order Lagrangian L(x,u, Vu), Noether's identity reads 



prMv(Z) + LDiC = Ei(L)(^ - v)?) + A 



r>T 



(3) 



where Ej is the Euler operator or variational derivative with respect to u 1 , |19j . Once stated, 
the verification of the identity is a straightforward computation. In the following sections, we 
will investigate how to use Noether's identity in the framework of elasticity, and apply the 
corresponding integral identities to establish non-existence results. The proofs are sketched, 
while the full details are left to the interested reader as exercises. 



3 Elastostatics 

We recall that the equilibrium equations for a homogeneous isotropic linearly elastic medium in 
the absence of body forces arise from the variational principle with Lagrangian 

1 1 1 n 1 ( n 

L Q (x,u,S/u) = -/i|| Vu|| 2 + -Gu + A)(V-u) 2 = -f,J2 (*4) 2 + 2^ + A ) I J2 u i 

i,j=l \i=l 

where the parameters A and /i are the Lame moduli. The squared norm of the deformation 
gradient matrix Vn refers to the sum of the squares of its entries, while V • u denotes the 
divergence of the deformation. The corresponding Euler-Lagrange equations are known as 
Navier's equations: 

LiAu+ (n + A)V(V -u) = 0, 

where the Laplacian A acts component-wise on u. Henceforth, we assume that /i > and 
fi+X > 0, thereby ensuring strong ellipticity and positive definiteness of the underlying elasticity 
tensor, [lOlfTT]. 

In this paper, we shall study boundary value problems for elastic bodies that are subject to 
a nonlinear body-force potential F(u). Thus, we modify the preceding Lagrangian 

L(x,u,Vu) = ^H|V-u|| 2 + \{» + A)(V • u) 2 - F(u), 

where we assume, without loss of generality, that F(0) = 0. The associated equilibrium Euler- 
Lagrange equations are 

fiAu + (n + A)V(V • u) + f{u) = 0, (4) 

where fi(u) = dF/du 1 are the components of the gradient of the body- force potential with 
respect to the dependent variables u. 

More generally, we consider Lagrangians of the form: 



L = -C%44-F(u), 



(5) 
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where again F(0) = 0, and 

e = - (Vu + Vn r ), with components e\ = -^{ u% k + u i)i 

is the strain tensor. The quadratic components in the Lagrangian ([5]) model the stored energy 
of a general anisotropic linearly elastic medium, while F(u) represents a nonlinear body-force 
potential. The elastic moduli C^J are assumed to be constant, satisfying 

= C% = C# = <%. (6) 

Thus in planar elasticity there are 6 independent elastic moduli, while in three dimensions 21 
independent moduli are required in general, [10] . Additional symmetry restrictions stemming 
from the constitutive properties of the elastic material may place additional constraints on the 
moduli. We may also assume 

Cga\a> > (7) 

for any matrix A = (dq). The less restrictive Legendre-Hadamard condition is that 

Ci-v^yj-WkWi > 

for any rank one matrix A = v <g> w. The Euler-Lagrange equations associated with ([5]) read 

Cg«i, + / i («) = 0. (8) 

In general, the most basic Pohozhaev-type identity is based on the associated Noether identity 
for the infinitesimal generator of an adroitly chosen scaling transformation group, [19] . 



Theorem 1. Let Q be a bounded domain in W 1 . Then the classical solutions of §8§ - that is 
u £ C 2 (fi) n C 1 (H) - subject to homogeneous Dirichlet boundary conditions on dfl satisfy the 
following Pohozhaev-type identity: 



^-^u k f k ( u )-nF(u) 



dx = -- [ Cfjvjya 3 l {x,u)ds, (9) 
1 Jan 



where v is the outward unit normal to dSX and (•, •) is the Euclidean scalar product in 
Proof. We consider the one-parameter group of dilations 

(x,u) > (Xx,X (2 - n)/2 u) 
with infinitesimal generator 
d 2 — n s d 



Xl dxi 2 du % 
According to JT|), ((2|), the first order prolongation of this vector field is 
cii d 2 — n • d n j d 

„ r V ; v — x 1 u yJ 

v dxi 2 du l 2 1 Qu\ 

Then one easily sees that 

pr( 1 ) V (L) + LDiC = ^ u k f k ( u ) - nF(u). (10) 

The identity ([9]) now follows from the divergence theorem using (J3j) , ©, (|10p . our assumption 
-F(O) = 0, and the homogeneous Dirichlet boundary conditions, taking into account that, on dfl, 

u{vi = u\v s . (11) 

See [22], p. 683] for more details on the last point. ■ 
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For the sake of completeness, we specialize the general elastic Pohozhaev identity to the 
isotropic case of the forced Navier equations in Q: 



n 



-u k f k (u) - nF(u) 



dx 



1 

2 .ion 



^H|Vu|| 2 + i(/J + A)(V-u) 2 



(x, u)ds, 



again subject to homogeneous Dirichlet boundary conditions on dQ. 

As a corollary, we obtain the following non-existence result. Recall that the domain Q is 
star-shaped with respect to the origin if (x, v) > for any x £ d£l. 

Theorem 2. Suppose that Q, is a star-shaped domain. Let the function 

F = F(s) = F(s l ,...,s n )eC 1 (R n ) 
satisfy the conditions F(0) = and 



n 



2 h dF 



s K O^-nF(s)>0, 



1. 



(12) 



for any s £ W 1 . We also suppose that the equality in (|12p holds if and only if s = 0. Then 
there is no non-trivial classical solution of the potential systems ([I]), ([8|), subject to homogeneous 
Dirichlet boundary conditions. 

Proof. This theorem follows easily from the identity Q, taking into account the positivity 
requirement ([7]) and star-shapedness condition. Indeed, any classical solution of ([8]) subject to 
homogeneous Dirichlet boundary conditions on d£l must satisfy the identity © . For ([7]) and the 
star-shapedness condition (x, u) > for any x € it follows that the right-hand side of (|9|) is 
non-positive. On the other hand, by (|12j) the left-hand side of ([9]) is positive unless u = in fL 
Hence u = 0. ■ 



4 Elastodynamics 

In this section, we turn our attention to hyperbolic elastodynamic systems of potential type: 



ul t + CMu{ l + f i (u) = 



y ij "-kl 



(13) 



in R x (1 with homogeneous Dirichlet boundary conditions on M x dQ. The corresponding 
Lagrangian is given by 



L = \C\]AA - \<< - F(u) = \c%uU - \a\u\ - F(u), 

where the second expression follows from the requirements ([6]) on the elastic moduli. 
Theorem 3. The classical solutions of the problem (]13p satisfy the following identity 

n — 1 

'U _ 



d 

dt ./ 



tE{u) + u\u\x k + n - l u i u i 1 



+ 



2 



dx 



-u k f k (u)-(n + l)F(u) 



dx 



ds, 



where 



E{u) = \ (C^eiel + «*«*) - F(u) = \ (Cgu\v> + 



F(u) 



(14) 



(15) 



is the energy density. 
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Proof. We introduce a vector field v which is the infinitesimal generator of the dilation group 

(t, x ,u) i— > (\t,\ x ,\( i - n y 2 u). 

The first order prolongation of v is given by 

cis d d 1 — rijd n + 1 • d n + 1 ■ d 

pr v = 'at + + 2 ^ av " — " — 

As a result, 

pr«v(L) + LAf = ^i fc / fc (u) - (n + l)F(u), (16) 

where the Lagrangian L is given by (|14p . Then, after some algebraic manipulations, the identi- 
ty (fT5|) follows from the Noether identity (|3|) combined with (JHJ), (fTTj) . (fT6|h the homogeneous 
Dirichlet boundary conditions, and, finally, the divergence theorem. ■ 



Let VL C R n be a ball of radius R centered at the origin. If we assume that u(t, x) decays 
sufficiently rapidly as R = \x\ — > oo, then the following conformal identity holds for the nonlinear 
hyperbolic system (fT3j) in R x R n : 



Corollary 1. The classical solutions of the problem (|13|) in R x R n that decay rapidly at large 
distances satisfy the identity 



di 



tE(u) + u\u\x k + U u l ul 



dx 



n — 1 



u k f k (u)-(n + l)F(u) 



dx. 



We observe that this result generalizes Morawetz's dilational identity for nonlinear wave 
equations, [HI [291 [30] , to elastodynamical systems. 

Finally, we consider a nonlinear hyperbolic system of so-called Hamiltonian type, [3], 



-u\ t + CVu> kl + H v 
-4 + Cflvja + H u 



0. 
0. 



(17) 



in R x with homogeneous Dirichlet boundary conditions on R x d£l. (The independent variab- 
le x must belong to an even dimensional space R 2m .) For such systems, we obtain a generalization 
of Morawetz's conformal identity |30j. 



Theorem 4. The classical solutions of the problem (|17p satisfy the following identity 

dx 



d_ 

di 



tE(u,v) + (x k u{v{ + x k v[u{) + n - 1 (au j v{ + bv j u{) 



n — 1 



[au 



dx 



+ 



_ H uk + bv k H vk ) 
{C§uivf + u\v\){x, v) + tC k l(u\v\v k + t^'u^) 



(is, 



(18) 



where the constants a and b are such that a + b = 2 and 



E(u, v) = C k Juivi + u\v\ - H(u, v). 
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Proof. In order to prove Theorem HI we use the same scheme as in the preceding Theorem [3j 
Namely, we consider a vector field v which is the infinitesimal generator of the dilation group 

(t,x,u,v) .— ► {Xt,\x,\ a ^- n ^ 2 u,X b{l ~ n)/2 v), 

where the constants a and b satisfy a + b = 2. Applying the first order prolongation 



pr«v 



+ 

to the Lagrangian 
1 



8 d q(l -n) i d b{l - n) i d 

t— + Zj- 1 U TT- H ~ V TTT + 



8t dXi 

a(l — n) 



a(l-n) \ t d 

—2 1 Ut M 







U 3 du l j 



6(1 -n) 



d_ 

dvl 



1 I v\— + 



2 ) idv) 



u\v\ - H(u,v) 



yields 



pr ( ^ 1) v(L) + LDiC 



n — 1 



, au 



k H u k + bv k H v k)., 



(19) 



(20) 



when a + 6 = 2. Then, after some additional work, the identity (|18p follows from (|19|) . (|2Up . 
(|lip . the homogeneous Dirichlet boundary conditions, and the divergence theorem. ■ 



Corollary 2. Let a, b and E be as in Theorem |H Then, provided u and v decay sufficiently 
rapidly at large distances, 



d_ 

dt 



tE(u,v) + (x k u 3 k v 3 t + x k v{u J t ) + 
n — 1 



n — 1 



av?v\ + bv-'ufy 



dx 



[au K H, l 



bv k H vk ) 



dx. 



Applications of these identities to the stability and scattering of waves in elastic media will 
be developed elsewhere. 



5 Further directions 

We emphasize that, in order to obtain Pohozhaev and Morawetz-type identities in elastostatics 
and elastodynamics by the Noetherian approach developed in [313], we have focussed our atten- 
tion on dilations, which are particular cases of conformal transformations. Further variational 
identities associated with other variational symmetries remain to be investigated. In particular, 
it would be interesting to analyze the variational identity for the semilinear Navier equations 
that corresponds to the first order generalized symmetry 

v= [vu) + (2» + \)d)4]-£j 

found in [18J. In fact, the variational and (at least in three dimensions) non- variational symme- 
tries for isotropic linear elastostatics were completely classified in [17\ [T8] and the systems not 
only admit point symmetries, but also a number of first order generalized symmetries. In the 
two-dimensional case, complex variable methods, as in [15] . are used to produce infinite families 
of symmetries and conservation laws. Also in the two-dimensional case, additional symmetries 
appear when 3^ + A = 0. In the three-dimensional case, when l\i + 3A = 0, Navier's equations 
admit a full conformal symmetry group, along with additional conformal-like generalized sym- 
metries. Although these restrictions are non-physical, they still lead to interesting divergence 
identities in the more general isotropic case, which can be applied to the analysis of eigenvalue 
problems, and also, potentially, the nonlinearly forced case. This remains to be investigated 
thoroughly. 
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